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Abstract 

For axisymmetric evolution of isolated systems, we prove that there 
exists a gauge such that the total mass can be written as a positive 
definite integral on the spacelike hypersurfaces of the foliation and the 
integral is constant along the evolution. The conserved mass integral 
controls the square of the extrinsic curvature and the square of first 
derivatives of the intrinsic metric. We also discuss applications of this 
result for the global existence problem in axial symmetry. 

1 Introduction 

The main unsolved problem in General Relativity is the global existence 
of solutions of Einstein's equations that describe the dynamics of strong 
gravitational fields. Symmetries have traditionally play an important role in 
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this problem. The presence of a symmetry reduces the degrees of freedom 
of Einstein equations and hence simphfy considerable its analysis. This is of 
course useful as a preliminary step to understand the full problem but also 
many models with symmetries have direct physical applications. 

In vacuum, due to Birkhoff's theorem, spherical symmetry has no dy- 
namics. For isolated systems, the next possible model with symmetries are 
axially symmetric spacetimes. It has been proved in [5j that no additional 
symmetry can be imposed to the spacetime if we want to keep the gravita- 
tional radiation and a complete null infinity. This result single out axially 
symmetric spacetimes as the only models for isolated, dynamical, system 
with symmetries. 

There exists many relevant physical models one can study in axial sym- 
metry: head-on collisions of two black holes, rotating starts and black holes, 
critical collapse of gravitational waves. These models have been studied nu- 
merically. Particularly relevant for the results presented in this article are 
the following references in which axial symmetry is imposed explicitly on the 
equations using cylindrical coordinates: [23], [2], [12], [H], [2^, pH], [30], [2H]- 

From the analytical point of view, there exists up to now no results for 
axially symmetric isolated systems (see the review [25] for results with other 
kind of symmetries in cosmologies). One of the purpose of this article is 
to initiate the study of this problem. We mention a related problem: cos- 
mologies with f/(l) symmetries. This symmetry was analyzed in [12], [H]- 
The equations are locally the same as in the axially symmetric case but the 
boundary conditions are radically different. In axial symmetry the Killing 
vector vanishes at the axis and that is the main source of difficulties. 

For axially symmetric data, the total ADM mass [Ij can be written as a 
positive definite volume integral over one spacelike hypersurface (see [H] , [S] , 
[22], [II])- This fact is likely to play a major role in the initial value 
(see the discussion in In order to study the implication of this integral 

formula for the evolution, the first natural step is to prove that there exists a 
gauge such that the mass integral holds not only at one hypersurface but in 
a whole foliation and it is conserved along the evolution. This is the subject 
of the present article. To understand this result, let us review the notion of 
mass in General Relativity and its conservation. 

The total mass of an isolated system is a boundary integral at infinity 
calculated on a given spacelike hypersurface. If we consider not only one hy- 
persurface but a foliation on the spacetime we can ask the question whether 
the mass is conserved along this particular foliation. That is, whether the 
boundary integral gives the same result if it calculated on different slices 
of the foliation. Any foliation is determined by the choice of a lapse func- 
tion and a shift vector. If the lapse and shift of the foliation satisfy some 
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fall-off conditions (this class is called asymptotically fiat gauges), then the 
mass boundary integral is conserved. This notions appears naturally in the 
Hamiltonian formulation of General Relativity (see [21] [1] [31] [32] ) . 

In any physical theory conserved quantities (in particular, conserved en- 
ergies) are very important to control the evolution of the system. However, 
in General Relativity, the conserved mass appears as a boundary integral and 
not as a volume integral (as, for example, in the wave equation). Hence it 
is not possible to relate the mass with any norm of the fields to control the 
evolution of them (for the wave equation the energy is precisely the norm 
of the wave). Axially symmetric systems represent a remarlsable exception. 
As it was mentioned above in this case it is possible to write the mass as a 
positive definite volume integral. However, this integral formula is valid only 
in a particular gauge. Moreover, this gauge can not be given a priori, it is a 
dynamical gauge which is prescribed as a solution of a system of differential 
equations. The natural question is if this particular gauge satisfy the fall-off 
conditions which guarantee the conservation of the mass. Our main result is 
that the answer to that question is yes. In other words, there exist gauge for 
which the mass can be written as a positive definite integral on each slice of 
the foliation and this integral is conserved along the evolution on this partic- 
ular foliation. Then, in this gauge, the conserved mass controls the norm of 
the metric and the second fundamental form along the evolution and hence 
this particular gauge is likely to be the most relevant one to study axially 
symmetric isolated systems (we further discuss this point in section [H]). 

The plan of the article is the following. In section [2] we describe our main 
result which is given by theorem 12. 1[ In section [3] we review the well known 
(2+l)-|-l formalism and compute the behavior of the fields at infinity and 
near the axis. In section H] we describe the gauge. We also derive the mass 
integral formula in the (2+l)-|-l formalism. This derivation is remarkably 
simpler than the ones presented so far in the literature. In section [5] we 
prove the main result. Finally, in section [6] we discuss the implication of this 
result for the evolution problem. 

2 Main result 

An initial data set for the Einstein vacuum equations is given by a triple 
{S, 'Jab, Kab) where S* is a connected 3-dimensional manifold, 7afc a (positive 
definite) Riemannian metric, and Kab a symmetric tensor field on S*, such 



3 



that the vacuum constraint equations 



DbK'''' - D^K = 0, (1) 
'^^^R-KabK^^ + K^ = 0, (2) 

are satisfied on 5*. Where D and ^^^R are the Levi-Civita connection and 
the Ricci scalar associated with ■jab, and K = Kab'j"'^- In these equations the 
indices are moved with the metric 'jat and its inverse 7"^. For simphcity, in 
the present article we assume S = M^. We expect, however, that the results 
presented here generalize to asymptotically euclidean manifold with many 
ends. 

We will further assume that the data are axially symmetric, which means 
that there exists a Killing vector field 7]°', i.e; 

£r,lab = 0, (3) 

where £ denotes the Lie derivative, which has complete periodic orbits and 
such that 

£,Kat = 0. (4) 

The Killing vector field has associated the following scalars with respect to 
the metric 'jat- The norm 

A' = 7abV^v\ (5) 

and the twist 

u' = eabcV^D'v'- (6) 

The prime in the notation of u' is justified in section 121 where we will show 
that uj' is related to the time derivative of the four dimensional twist potential 

UJ. 

The data are called asymptotically flat if there exists a coordinate system 
x'^ (/i = 1,2,3) in M.^ such that the metric components in these coordinates 
satisfy the fall-off conditions 

7f,u = + Ooo(r"^/^), (7) 

and the extrinsic curvature components 

K^. = Ooo(r-3/'), (8) 

where we write / = Ooo{r'') if / satisfies (9"/ = o(r'^~l"l), for all a. Here 
6^1, is the euclidean metric in the coordinates x^, r is the euclidean radius 
r = x^x"b^y, d denotes partial derivatives with respect to the coordinates 
and a is a mult index. 
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The data are called maximal if the trace of Kab is equal to zero, i.e. 



K = r'K^b = 0. (9) 
Following [IH], we define the two dimensional quotient metric qab by 

Qab = lab (10) 

In [H] it is proved that for any axially symmetric and asymptotically flat met- 
ric 7afe there exists a coordinate system (p, z, 0), which we call the isothermal 
coordinates, such that the metric Qab and the Killing vector have the following 
form ^ 

q = e''^idp^ + dz^), Tr= , (11) 



and the norm is given by 

A = pe^l^ (12) 

where u and a are smooth functions. To these coordinates we associated 
Cartesian coordinates {x,y,z) via the standard formula x = pcos0, y = 
psin0. It follows from the results proved in [ll| that in the coordinates 
{x, y, z) the metric 7afc is asymptotically flat. Hence, without lose of general- 
ity, in the following we take to be (x, y, z). 

In isothermal coordinates we consider the following integral 
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m = — / dz 
16 



p2u 12 

\da\' + ^-^ + 2e'-K'^'K,b 



pdp, (13) 



where \da\ = afp + a^^. Note that the volume element pdzdp is equivalent to 
the volume element in Cartesian coordinates dxdydz where the integration 
domain is M^. We will see in section H] that for asymptotically fiat, axially 
symmetric and maximal data m defined by (fT3ll is the total mass. 

The mass integral (JT3l) is positive definite, and is zero only for fiat space- 
time. The first two terms in the integrand contain the square of first deriva- 
tives of the metric and the third term contains the square of the extrinsic 
curvature. 

We want to prescribe gauge conditions for the evolution such that the 
integral (|T3i) holds not only at the initial surface but in any surface of the 
foliation. We prescribe the lapse by the requirement that the maximal con- 
dition is preserved under evolution. The shift vector is prescribed by the 
condition that the isothermal coordinates are preserved by the evolution. In 
this way we obtain a coordinate system {t,p,z,(j)), we call it the maximal 
isothermal gauge. These gauge conditions are explicitly given in section |H 
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This gauge has been used in numerical evolutions in the following references 
[23], [12!, [S], [2n], [2H]- This is a dynamical gauge, namely, the lapse and 
shift are calculated in terms of the initial data at each step of the evolution. 
The main result of this article is given by the following theorem. 

Theorem 2.1. Let (M'^, 7^6, i^a^, r/'*) be an axially symmetric, maximal, vac- 
uum, asymptotically flat initial data. Then, there exist a unique solution of 
the maximal-isothermal gauge equations and the mass integral given by ( |T3l) 
is conserved along the evolution, that is 

^ = 0. (14) 
dt ^ ' 

This theorem is proved is section HI We divide the proof in lemma 15.11 
theorem 15.21 and theorem 15.31 For the sake of clarity, in this section we 
presented the results using the standard 3+1 decomposition. However, for 
axially symmetric spacetimes there exists a more natural decomposition, the 
(2+l)+l formalism described in the section[31 The maximal isothermal gauge 
conditions and the mass formula are better expressed in this formalism. 



3 Axisymmetric evolution equations 

In this section we review the (2+l)+l formalism, which has two parts. First, 
a reduction of the field equations by the action of the symmetry [18] and 
second a time plus space decomposition of the reduced equations (see [23] . 
[H] , [26] ) . We also give some useful formulas which relate the quantities in the 
(2+l)+l formalism with the natural ones in the standard 3+1 decomposition. 



3.1 The (2+l) + l formalism 

Consider a vacuum solution of Einstein's equations, i.e., a four dimensional 
manifold V with metric gab for which the Ricci tensor ^^^TZat vanishes. Sup- 
pose, in addition, that there exists a spacetime Killing vector 77°. We define 
the norm and the twist of 77°, respectively, by 

A' = V'^v'gab, UJa = eafccdr/'VV, (15) 

where Vq is the connection and eatcd the volume element with respect to gab- 
Assuming that the manifold is simply connected and using ^^^TZab = it is 
possible to prove that Ua is the gradient of a scalar field uj 

Ua = VaUJ. (16) 
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In our case the Killing field will be spacelike, i.e. A > 0. 

In the presence of a Killing field, there exists a well known procedure to 
reduce the field equations [H]. Let M denote the collection of all trajectories 
of ?7", and assume that it is a differential 3-manifold. We define the Lorentzian 
metric hab on M by 

Qab = flab + (17) 

Einstein vacuum equation are equivalent to the following set of equations 
intrinsic to (AA, hab) 

V.V^A = -^V"a;V,^, (18) 
V„V"a; = jW^'ujWaK (19) 

^^^Tlab = jVaVbX + ^ (VaOuVbUJ - VcOoV'uhab) ■ (20) 

where Va and ^^^TZab are the connexion and the Ricci tensor of hab- The 
system of equations (fT8l) - (!20l) can be interpreted as Einstein equation in 3- 
dimensions (equation ( l20i) ) coupled with a matter sources given by equations 

We make a 2 + 1 decomposition of (A/", hab) and the field equation fl20l) . 
Let n°' be the unit normal vector orthogonal to a spacelike, 2-dimensional 
slice S. By construction we have n°'ria = 0. And hence the norm defined in 
( |T5l) is equivalent to ([5]). The intrinsic metric on S is denoted by qab and it 
is given by 

hab = -TiaUb + Qab- (21) 

Our convention for the signature of hab is ( — I" +)• In terms of the lapse a 
and shift vector the line element takes the form 

h = -a^dt^ + qij{dx' + I3'dt){dx^ + 13^ dt), (22) 

where i, j = 1, 2. The extrinsic curvature Xab of the slices S is given by 

Xab = -ql^crib = ~^£nqab- (23) 

We write the field equations ( j20l) as evolution equations for {qab, Xab) 

Qab = -2aXab + £(}qab, (24) 
Xab = £(3Xab + a {xXab + ^"^^ Rab "^^^ T^ab - '^■XacXt) " DaDbOi, (25) 



^opposite sign convention with respect to 
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and constraint equations 

^'^i?-x'^'Xa6 + x'=/i, (26) 
D'^Xab - DhX = Jb. (27) 

Here a dot denotes partial derivative with respect to the t coordinate, ^'^^Rat 
is the Ricci tensor and Da is the connexion with respect to Qab, ^^^T^ab denotes 
the projection of ^^^TZab into 5* and 

fi = 2^^'^nabrfn^ +(^) 7^, (28) 
J, = -qln^^^n^. (29) 

The scalar /x, which appears as 'matter source' in the Hamiltonian constraint 
fpS|) . will be relevant in the computation of the mass. We will need its explicit 
expression 

^ = 2^ + 2Da [-^j + + ^ K + l^^l ) ' (30) 

where the prime denotes directional derivative with respect to n", that is 

A' = n'^VaA = £n\ (31) 

Note that 

A' = a-\\ - £p\). (32) 

To obtain equation (IHUI) we have used equations (fT^ - (fT^ . 

Finally, it convenient to decompose the extrinsic curvature Xab in its trace 
X and trace free part kab, that is 

kab = Xab - Qab^- (33) 

3.2 Relations with the 3+1 decomposition 

Take the definition of the spacetime twist potential (1151) and (1161) . The 
derivative u' is given by 

a;' = n"V„^ = n'^eabcdri'V^'r]''. (34) 

Note that rf-eabcd is the volume element on the 3-dimensional slice S de- 
noted by eabc in section [21 The covariant derivative V a in ([MD appears 
antisymetrized and hence we can replace it by any covariant derivative, in 
particular Da- Then, we obtain that (IMl) is equivalent to ([6]). 



8 



The twist potential u can be calculated in terms of the initial data as 
follows. Define the vector S"' by 

Sa = KabV'- (35) 

Using the definition of Kab 

Kab = -jycUb, (36) 

we obtain 

Sa = -n'V.Va, (37) 
where we have used the Killing equation V(aVb) = 0, equation rjau"- = and 
the following expression for the metric jab 

gab = lab - UaUb. (38) 

We use the following expression for the covariant derivative of rja (see 

m) ^ ^ 

^aVb = ^^^abcdV''^'^ + -^VlaDbjX- (39) 

Using (139|) and (137|) we obtain 

eafcc5V=-^^a^- (40) 

Note that the left hand side of equation (HOi) is calculated only with the initial 
data. This is the desired equation. 

The relation between the trace of the extrinsic curvature K of the slice 
S and the trace of extrinsic curvature x of the slice S is given by 



If we prescribe K = then we have 



K = X-j. (41) 



X=j. (42) 



A 

The relation between the two extrinsic curvatures is given by 

Xab = q'LqtKcd. (43) 

The following component of the extrinsic curvature will be used in the next 
section 

K^^ = KabV^v' = -A'A. (44) 
Using equations ( |35l) . (HOl) . (H3ll . and (jSj) we obtain the following expression 
for the square of Knh 



K^'Kab = X^'Xab + ^ + (45) 
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4 Maximal- isothermal gauge and mass inte- 
gral 

In this section we review the well known maximal-isothermal gauge (see [23] , 
[T7j [H], [2S], [2H]) using the (2+1 )+l formalism presented in the previous 
section. We analyze the fall off behavior of the fields in this gauge near the 
axis an at infinity. Finally, we present a new derivation of the mass integral 
formula. 

4.1 Gauge 

The maximal gauge condition is well known. In the standard 3 + 1 decom- 
position, the equation for the lapse is given by 

A^a = aK^^Kab, (46) 

where denotes the Laplacian with respect to the metric 7af,. 

The shift vector is fixed by the requirement that the isothermal coor- 
dinates are preserved under the evolution. In isothermal coordinates, the 
2-dimensional metric qab has the form (|TT1) . If we take a time derivative to 
ffTTj) we obtain 

Qab = 2U Qab, (47) 

that is, the tensor g^fe is pure trace. Hence, the trace free part (with respect 
to Qab) of the evolution equation is given by 

{Cg(3)ab = 2aKb, (48) 
where Cq is the conformal Killing operator with respect to Qab 

iCg/3)ab = Da/3b + Db/3a - QabD^/S''. (49) 

Equations ( HHl) constitute a first order elliptic system which determine (3°' 
under appropriate boundary conditions, as we will see in section [51 Summa- 
rizing, the lapse and shift for the maximal isothermal gauge are determined 
by equations fj^ and fj^Hj) . 

4.2 Fall off behavior and axial regularity in isothermal 
coordinates 

The fall of conditions ([7]) for the 3-dimensional metric implies that (see [E]) 

u = a = Ooo(r-^/'). (50) 
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From equations fH5l) and equation ([H]) we obtain 

X^, = Ooo(r-=^/'). (51) 

The function A'/ A will be important in the next section, we need to compute 
its fall off behavior. This behavior is not an obvious consequence of ([H]) 
because the relation between these two functions involves a quotient by p^. 
To compute the fall off of this function we proceed as follows. In Cartesian 
components we have 

r/^ = -y, = x, (52) 

and hence 

K^^ = y^K^^ - 2xyK^y + x^Kyy. (53) 

Following [Hj it is enough to consider the plane x = which is transversal to 
the Killing vector. If we use ( l44l) . ( l53l) and ( |T2l) and then evaluate at x = 
we have 

- y = ^ = e-^K,^. (54) 

Hence, at this plane the function A'/A is smooth and it has the same fall off 
as Kah, namely 

y = o^{t-'''). (55) 

Since the plane x = is transversal to the axial Killing vector it follows that 
fl55p holds everywhere. 

The regularity of the metric 7afe and the extrinsic curvature Kab at the axis 
implies restrictions on the behavior of the different quantities (for a detailed 
discussion of this issue see [26j and [29j). We summarize here two important 
consequences of axial regularity. First, the function q defined by 

g = M-|, (56) 

vanished at the axis (see [T^) 

q{p = 0) = 0. (57) 

Second, the functions involved satisfy parity conditions with respect to the 
p coordinate. In particular we have 

u,q,a,kp are even functions of p, (58) 

and 

A, are odd functions of p. (59) 
Note that odd functions vanishes at the axis. 
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4.3 Mass formula 

Consider the Hamiltonian constraint ( !26l) . Using equations ( fT2l) and ( !56l) we 
write the Hamiltonian constraint in terms of p, a and g in such a way that 
there are no singular terms at the axis. To do this, we use the following 
elementary identities. 

The Ricci scalar of qah is given by 

where A is the flat 2-dimensional Laplacian 

Au = u^pp + u^zz- (61) 
Using equation (fT2l) we obtain 

--\Dlog\\^=\D{^ + \ogp)\''. (62) 



A2 

Finally, the p coordinate satisfy 

(^)Alogp = 0. (63) 
where we have defined ^^-"A as the 3-dimensional flat Laplacian 



(3)A« = Au + ^. (64) 
P 



Using equations ( l56l) . ( 1601) . ( l62l) and ( l63i) the Hamiltonian constraint ([26] 
can be written in the following form 



^2u 



2x''\ab + 2x(-X + 2y) + \Da\^ + ^ (a;'^ + \Du\^) 



(65) 

We define the mass as the integral of the right hand side of this equation, 
namely 

m = ^ ^ (2k'^'K, -X' + 4Xy + \Da\' + ^ {u" + iZ^c.^) j prfV„ (66) 

where we have used (l33l) and rfV^ = e^^dpdz denote the volume element with 
respect to qah- Note that we introduce a weight p in the integral. If we 
impose the condition K = Q (using equation (HTI) ) we get that the integrand 
is positive definite, namely 

m = ^ ^ (^^"'^"^ + + '^"^'^ + ^"^'^ + '^"^''0 ^^'^^ 
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For maximal 2- dimensional slices x = the integral flbbl) is also positive 
definite. However, this choice will not give an equation for the lapse function 
(analog to ( H6l) ) with positive definite solutions. To prove the equivalence of 
equation ( 1671) with ( |T3|) we use ( H5l) . 

The crucial point is that this volume integral can be written as a boundary 
integral at infinity using the left hand side of equation We have 



m 



/'OO /"OO 

- / dz {^^^Aa + Aq) pdp. {6t 

4 J-oo Jo 



This integral can be converted into a boundary integral at infinity as follows. 
For the first term we use that pdpdz is the volume element in M? and then 
we can use the divergence theorem in three dimension to obtain 

/oo 1*00 
dz / ^^^Aa pdpdz = lim / drcrrpdO. (69) 
-oo Jo Jo 

The second term in (lUHj) can be also written in divergence form 

j dpdz {q^pp + q^zz)p = j dpdz {{pq^p - q)^p + {pq,z),z) ■ (70) 

We use the divergence theorem in two dimensions to transform this volume 
integral in a boundary integral. Namely, let VL be an arbitrary domain, we 
have 

dpdz {{pq^p- q)^p + {pq^^)^^) = (t V-nds, (71) 
Jn Jan 

where n is the 2-dimensional unit normal, ds the line element of the 1- 
dimensional boundary dQ and V is the 2-dimensional vector given in coor- 
dinates (p, z) by 

V=iipq,p-q),ipq,z)). (72) 

Let Q be the half plane p > 0. By ( 1571) and the assumption that q is smooth 
we have that the vector V vanishes at the axis. Then the only contribution 
of the boundary integral is at infinity, namely 



A2qpdpdz = lim / (rdrq — q) pd9. (73) 
p>o Jo 

Summing (I7H|) and (lU^ we obtain the final expression for the mass as bound- 
ary integral 

m = — lim / (dr(T + drq—-)rpd6. (74) 

4 r^oo Jq \ r / 
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We can write this expression in terms of u instead of q 

4 r-»oo Jq y2 2r r J 

Up to now we have not mention the relation of m with the total mass. In fact 
in order to prove theorem 12.11 we do not need any other information. The 
quantity m is a positive definite integral that can be expressed as a boundary 
integral at infinity. As a consequence we will prove that it is conserved along 
the evolution in the maximal-isothermal gauge. 

However, it is possible to show that m is precisely the total mass of the 
spacetime. There is an interesting and subtle point here. The equivalence of 
this formula with the total mass have been previously proved using a strong 
fall off condition for q (see [Oj, ^19j, [IS]), namely 

g = Ooo(r-^/'), (76) 

which implies that q does not appears in the boundary integral (1731) . Re- 
cently, this equivalence was proved without this restriction [H]. In the de- 
duction we have made, we have not use any restriction for q which is consis- 
tent with [Hj. Note the simplicity of this deduction (which is much in the 
spirit of the energy defined in P^2] for cosmologies) compared with the ones 
using three dimensional slices like [19] . 

We also mention that in the next section we will prove that the strong 
decay condition ( 1761) is preserved by the evolution. 



5 Existence of the gauge and mass conserva- 
tion 

We begin with the lapse equation fj^6|) . This condition does not require 
any symmetry on the data. The existence and uniqueness of the solution 
of this equation in the standard 3 + 1 decomposition are well known. For 
completeness we review this result. 

Lemma 5.1. There exists a unique smooth strictly positive solution a of 
equation ( l46l) with the following fall off 

a-l = o^{r-"^). (77) 

Proof. Set 

a = a — 1, (78) 
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then, using fj^Uj) . we have that 

A,a -fa = f, (79) 

where / = K"'^Kab > 0. The fall off behavior assumed for Kah in ([H|) implies 
that / G Hl^i^ for all s. Here denote weighted Sobolev spaces. See [3] [7j 
[in] and reference therein for definition and properties of these spaces. We 
use the index notation of [S]. 

The decay assumed on the metric 7^;, ([7]) implies that the operator A^ — / : 
— *■ i?|_2 is an isomorphism for 5 = —1/2 (see for example pjLj, [21]). 
Note that we use this theorem for dimension 3. Then there exists a unique 
solution d G '^^ equation (ffOll . 

Since a — as r ^ 00 (because a G H^^^2)j it follows that a — 1 as 
r ^ c». Hence, we can use the version of the maximum principle for non- 
compact manifolds presented in pjj to conclude that there exists e > such 
that a > e. 

□ 

If the data is axially symmetric then the solution a will be also axially 
symmetric: take the Lie derivative with respect to rj"' to both sides of equation 
(^B, we get 

A^{£r,a) - f {£r,a) =0, (80) 

since G (^o see this, take the expression of 77 in Cartesian coordi- 

nates 0521) ). the isomorphism theorem mentioned above implies that £^q; = 0. 
Finally, we note that since a is smooth in {x,y,z), it follows that an axially 
symmetric a is an even function of p. 

We consider now the equation for the shift. Instead of solving the first 
order system (HHll we will solve a second order system obtained from taking 
a derivative to fHHl) . Then, we will show that a solution of the second order 
system is also a solution of the original equations (jUj) under appropriate 
boundary conditions. 

Under the conformal transformations (fTTjl , the conformal Killing operator 
rescale like 

jCgif3)ab = e2"£(e-2"/3),,, (81) 
where C is the flat conformal Killing operator 

{C(3),b = daPh + dh(3a - qabdcP'. (82) 
We take a flat divergence to equation fjlHl) and use equation (IRTl) to get 

A{e-^^Pa)=2d\ae-'^Kb), (83) 
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where we have used 

d\C(3U = ^Pa. (84) 

If we contract equation (l83l) with and use q""^ = e~^"(5"^ we obtain our 
final equation 

Ap" = 2d\akl), (85) 

where = q°'^l3i, and = q'^^kcb 

As it was mentioned in section 14.21 the smoothness of the metric at the 
axis imphes that the relevant functions satisfies parity conditions in the p 
coordinate. In particular, the components of the shift vector j3°' should satisfy 

/3^(p, z) = r (-P, z), P^ip, z) = -pP{-p, z). (86) 

Theorem 5.2. There exists a unique, smooth, solution (3°" of equation fl85p 
with the following fall off behavior 

= Ooo(r-i/2). (87) 

The solution satisfies the parity conditions ( l86i) and the gauge equation ( l48l) . 
Moreover, the quotient P^/ p is smooth and have the following fall off 

^ = Ooo(r-^/^). (88) 
P 

Remark: the fall-off condition (IHHjl will be essential in the proof of the 
mass conservation theorem 15.31 



Proof. We distinguish the two components of equation ( l85l) 

ApP = PP, Ap' = F\ (89) 

where 

PP = 2d^{akP), P' = 20^ {ak]). (90) 

The function a is even in p, using ( |58l) and (1591) we obtain that PP is odd 
and P^ is even in p. 

In terms of the cylindrical coordinates (p, z), the source functions P are 
defined only for p > 0. But because they satisfy parity conditions we can 
smoothly extend them to all by the following prescription 

PP{p, z) = -PP{-p, z), P^p, z) = P^{-p, z). (91) 

Hence, we can consider equations flH^ as Poisson equations 

Au = /, (92) 
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in with decay conditions at infinity flH7|) . Uniqueness of solution of 
under the decay condition flH7|) follows immediately integrating by parts the 
Laplace equation. Namely, for any domain Q we have 



uAu = / d{udu) — \du\ = (h udnU — / \du\ , (93) 
n Jo. Jan Jn 

and hence for a solution of An = we have 



udr,u= / \du\\ (94) 
an Jn 

Let f2 be a ball of radius R, and take the limit R —* oo. By the decay 
condition (!87|) . the boundary integral in vanishes. Then, it follows that 
M must be constant, by the decay condition it must be zero. Since the 
homogeneous equation has only the trivial solution uniqueness follows. Also, 
for solutions in this class, integrating equation ( |92l) we obtain that the source 
must satisfy the condition 

' / = 0. (95) 



We need to verify for our source functions. The function is odd in 
p, hence it satisfies ( lU^ automatically. For F^ we use that it can be written 
in divergence form 

Then we have 



F' = div\ v' = 2aki. (96) 



F' = 2 F' = 2 div' (97) 

Jp>0 Jp>0 

using the divergence theorem in 2-dimension we get 

div' = / vP\p=odp+ lim (b v^'ds, (98) 
p>o Jo ^-^"^ Jcr 

where Cr denotes the semicircle with r = R. By the decay condition fl^T]) 
for kab and ( !77|) for the lapse, we have that the second boundary integral 
vanishes. For the first one, we equation ( l59i) to conclude that 

v'\p=o = 0. (99) 



And hence we obtain the desired result that F^ satisfies ( 195|) . 

To actually obtain a solution of equation ( l92i) we use the Green formula 
for the Laplacian in two dimension G{x,x') = — log|x — x'p. A solution u 
has the integral representation 

u{x)= [ G{x,x')f{x')dx'. (100) 
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In polar coordinates (r, 9) (related to (p, z) by the standard formula z 
r cos 6 p = r sin 6) we have the well known expansion of G 



G{r,e]r',e') = -21ogr + 2^ — cos(n(^ - 6')) for r > r', (101) 

n=l 

°° n 

= -2 log r' + 2 ^ ^ cos(n(0 - 9')) for r < r'. (102) 

n=l 

We use the trigonometric identity 

cos(n(^ - 6')) = cos(n^) cos(n^') + sm{ne) sm{ne'), (103) 

to split the sum in terms which are odd functions of 6 (and hence of p) and 
terms which are even functions of 6. 

If the source function / in the integral f llOOl) is odd in p, then the terms 
involving cos{n6') vanish, for even source functions the terms with sin (n6'') 
vanish. 

We are interested in the fall off behavior of the solution. It is convenient 
to split the source functions F in two term, one with compact support in 
some ball and the other which pick up the decay behavior. Let x : M — > M be 
a cut off function such that x e C~(M), < x < 1, x(^) = 1 for < t < 1, 
X(t) = for 2 < t and write Xnir) = xir/R), /oo = (1 - Xr)/, where R is 
an arbitrary positive number. 

Using the expansion fllOll) we compute the integral fllOOp . Note that the 
term with log r vanish because ( l95l) . For r > R we have for 

= 2f2'-^^ [r-M„ + B^ir) + C„(r))] , (104) 



n 

n=l 



where 



rR 

An= r"'sm{n9')FPr'dr'd9', (105) 
Jo 

B4r) = r-'' r''''sm{ne')Fl^r'dr'de', (106) 
Jr 

r°° 1 

C^{r) = r'^ —sin{ne')F^r'dr'de'. (107) 

J r ^ 

For we have a similar expansion for r > i? 

(3^ = 2Y^ [r-M„ + i?„(r) + C„(r))] , (108) 



n 

n=l 
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where 

nR 

An= I r"'cos{n9')F'r'dr'de', (109) 
Jo 

S„(r)=r-"/ r"'cos{ne')F^r'dr'de', (110) 
Jr 

poo ^ 

C^{r)=r'' — cosine') F'r'dr'de'. (Ill) 
Jr r'"- 

We use the decay condition fl^ to obtain that the sources satisfies 

FP = F' = Ooo{r-^^^). (112) 

From (11121) we get that 

Bn{r) = Cn{r) = Ooo(r-i/2). (113) 

These are the solutions of our problem. It is clear that they satisfy the parity 
conditions (!86ll . 

The important property of the representation fll04p is that it allow us to 
prove flHHl) . In effect, we use the trigonometric identity [20j 



n-l 

sm(n9) = 

k=0 



■ J2 Q(cos^)'=(sin^)"-'=sin Q(n- A;)7r^ (114) 
to obtain that the function 9„ defined by 

Qnie) = = J2 (t) (cos^)'=(sin^)"-^-^in (^-{n - k)n] (115) 

is an smooth function. Hence, using (11041) . we obtain the expansion 

^ = 2 V ^ [r-"+iA„ + r-ii?„(r) + r-'C^{r))] . (116) 

' n=l 

Using a similar argument as above we obtain (|88l) . 

Finally, we need to verify that our solution of equation (l85l) is also a 
solution of the first order system (HHl) . Define the trace free tensor tab by 



U = {{CqP)ab-2akab)e-^''. (117) 

We have proved that 

dHab = 0. (118) 
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Let tpp = ti and tp^ = t2, then equations (11181) are given by 

dih + d2t2 = 0, dih - d2ti = 0. (119) 

These are the Cauchy-Riemann equation for a complex function / = ti + it2. 
That is, a non trivial solution of flll9p implies that / is an entire function 
on the complex plane. But ti and t2 decay to zero at infinity. Hence, by 
Liouville's theorem, / vanishes. And then, we obtain that (3 is also a solution 
of the gauge equation f HHl) . □ 



Solutions of equation fl92|) can be obtained by other methods, for example 
see Proposition 2.6 in |14j. In the proof above we have used the Green 
function in order to prove ( IHHi) . It is likely that also this kind of expansions 
can be obtained by other methods like the ones used in [T3] . 

To prove the mass conservation we will make use of the evolution equation 
given by the trace of ([211), namely 

2u = -ax + Dap". (120) 

In terms of partial derivatives this equation is written like 

2u = -ax + dip' + 2p'diU. (121) 

Theorem 5.3 (Mass conservation). We have 

^ = 0. (122) 
dt ^ ' 

Proof. We use equation fl7^ to get 



dt 



lim / {]-dra+ ^ + drU--] rpde. (123) 
16 r^oo 2r r J 



In order to prove that this integral is zero we have to compute the decay of 
the integrand at infinity. We begin with the terms containing a. 
We use the relation 

A' a' p' , , 

and 

p' = - — , a' = - {& - P'd.a) , (125) 
a a 

to conclude that 

A' BP 
a = 2a- + p'dia + 2—. (126) 
A p 
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Using equations (1331) and flHHj) we obtain 



a = o. 



'oo 




(127) 



For ii we use equation (I12ip and the decay on [3 to get 



u = o. 



'oo 




(128) 



Using 01271) and (11281) we obtain that the boundary integral fll23p vanishes. 



We mention also another consequence of equation fll2ip . From this equa- 
tion we deduce that 



Then, if the function q has initially the stronger decay (1761) . this decay will 
be preserved by the evolution. 

6 Final comments 

In this final section we discuss the implication of the results presented here 
for the evolution of axially symmetric isolated systems. 

The first question one need to face in the evolution problem is the choice 
of gauge. In axial symmetry there have been studied different kind of gauge 
conditions (see for example [2] [22] )• The mass conservation formula single 
out a particular one. If we define appropriate Sobolev norms, then the mass 
formula essentially implies that the norm of the metric and the norm 
of the extrinsic curvature are bounded along the evolution. This is very 
desirable property which is not present in other gauges. The mass conserva- 
tion formula strongly suggests that the most convenient gauge for the axially 
symmetric evolution problem is the maximal isothermal one. 

As it was mention in the introduction, the conservation of mass is closely 
related with the Hamiltonian formulation of General Relativity. The maximal 
isothermal gauge is a dynamical gauge that depends on the time coordinate. 
The analysis of time dependent gauge conditions and its relation with con- 
served quantities in the Hamiltonian formulation was recently studied in [32] 
|31j . The maximal isothermal gauge satisfies the fall off conditions described 
in these references. 

In axial symmetry, even when the gauge conditions are fixed there ex- 
ists many possibilities to extract from Einstein equations a set of evolution 
equations. The reason for this ambiguity is that Einstein equations in 3- 
dimensions has no dynamics, this essentially means that equation ( !20l) can 



□ 




(129) 
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be replaced by the constraint equations fl^ - flTTl) (see If we chose to 

do so, then we obtain an evolution scheme in which the evolution equations 
are given by (fT8l)-( fT9|) . and the other equations (including gauge conditions) 
are elliptic constraint. But there exists others alternatives. Following [26] , 
we can classify them by the number of evolutions equations used vs the 
number of elliptic constraint equations. The scheme presented above has 
the minimum number of evolution equations and the maximum of constraint 
equations. This is called a fully constrained system. On the opposite side we 
have a system in which we do not solve for the constrains and we solve the 
evolution equation (|2^ - ( |25l) . This is called a free evolution scheme. In be- 
tween we have other possibilities to construct partially constrained systems. 
In all cases, the gauge conditions are solved as elliptic constraints. That is 
all schemes have a mix of evolution and elliptic equations (see the discussion 
in chapter 3 of |26j). 

To prove the mass conservation we have used two equations: the Hamil- 
tonian constraint (l26l) and the evolution equation for the conformal factor of 
the 2-dimensional metric fll20p . In a free evolution scheme the Hamiltonian 
constraint is not solved. Hence, the mass will not be conserved for arbitrary 
data. It will be only conserved for data that satisfies the constraint equa- 
tions. Then, in this case the mass integral formula will be not useful for 
controlling the evolution. The same will happen with partially constrained 
systems in which the Hamiltonian constrain is not solved or the evolution 
equation fll20p is not used. Such systems was discussed in [2S] and in [H] 
(in this reference, this kind of systems appears when the evolution of the 
conformal factor is used instead of the Hamiltonian constraint). 

On the other hand, for a fully constrained scheme (like the ones used in 
[H], [2H]) and a partially constrained scheme in which both the Hamiltonian 
and the evolution equation for the metric are used (like the one studied in 
|17j ) the mass will be conserved. The mass conservation formula single out 
these two schemes. The natural problem now is to study their well possness. 

Finally, it is interesting to note that the integral mass formula has direct 
application in numerical simulations. In numerical simulations of isolated 
systems, one often uses a position-dependent resolution that is high in the 
central region and much lower close to the outer boundary of the compu- 
tational domain. Hence computing the mass as a volume integral is more 
accurate than as a boundary integral. The resulting approximation to the 
mass is much better conserved during the evolutioiu [27 \ . 

thank O. Rinnc for pointing this out to me 
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